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An additional Z(, symmetry hidden in the fermion and Higgs sectors of the Standard Model has 
been found recently. It has a singular nature and is connected to the centers of the SU (3) and SV (2) 
subgroups of the gauge group. A lattice regularization of the Standard Model was constructed that 
possesses this symmetry. In this paper we report our results on the numerical simulation of its 
Electroweak sector. 
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I. INTRODUCTION 

It is the conventional point of view that all the symmetries of the Standard Model (SM), which must be used 
when dealing with its discretization, are known. Recently it was found that there exists an additional = Z2® 
symmetry in the fermion and Higgs sectors of the SM l] . It has a singular nature and is connected to the centers Z3 
and Z2 of the SU{i) and SU{2) subgroupS|l^. The gauge sector of the SM (in its discretized form) was redefined in 
such a way that it has the same perturbation expansion as the original one, while keeping the mentioned symmetry. 
The resulting model differs from the conventional SM via its symmetry properties. Therefore, we expect it would 
describe nature better than the conventional discretized Standard Model, if the additional symmetry does takes place. 

It is worth mentioning that the present status of the Standard Model on the lattice implies that it must be considered 
as a finite cutoff theory This is in agreement with the understanding that the SM does not describe physics at 
extremely small distances. Hence, it is sufficient to consider a cutoff A that is finite but much larger than all observed 
energies. The consideration of the infinite cutoff limit would be an attempt to continue the Standard Model to 
infinitesimal distances. Now it is believed that this attempt leads to a trivial continuum theory ^J- Nevertheless, at 
energies much less than the cutoff we can calculate any physical variable. 

So we can examine our model in order to understand whether the considered additional symmetry is important 
for the discretization of the Standard Model or not. As a first step in this direction we investigate numerically the 
quenched Electroweak sector of the constructed discretized SM. 

II. A HIDDEN SYMMETRY 

In this section we repeat our construction reported in in continuum notation. This is done in order to demonstrate 
the universal (regularization - independent) nature of the additional symmetry found. 

A. The Standard Model 

The Standard Model contains the following variables: 

1. The gauge fields associated with the symmetry group SU{?>) x SU{2) x U{1), which are the elements of the 
corresponding algebras: 

= Z^Xaesu{3), 

Bieu{l) = (-oo,cx3). (1) 

(Here, Xa are the Gell-Mann matrices, and at are the Pauh matrices.) The corresponding SU{3), SU{2), and U{1) 
field strengths are: 

Hij = d[iZj] + i[Z^, Zj], 
Gij = d[iAj] +i[A.i,Aj], 

= (2) 
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2. Anticommuting spinor variables, representing leptons and quarks: 

e n T \ I u c t 



Vf. Vt, J ' \ d s b 



(3) 



3. A scalar doublet 

a = l,2. (4) 

The action has the form: 

S = Sg + SH + Sf, (5) 

where we denote the fermion part of the action hy Sf, the pure gauge part is denoted by Sg, and the scalar part of 
the action by Sh- 

As usual, we consider Sg in the form: 



Sg = - I d'^x 



355(7(3) 2ffs;7(2) 9u{l 



(6) 



where we introduced the gauge couplings gsu{i)^ 9su{2)) and gu{i)- 
The scalar part of the action is 



Sh= j d^x\{d^ + iA^+iB^)^'' + J d^xV{\^\), (7) 

where l^d^"]) is the potential, which has a minimum at a nonzero value of ^ = v, causing spontaneous symmetry 
breaking. 

We express Sf through left-handed doublets L and right-handed singlets R of fermions: 



I ^ 1-75 / T 



3 - 2 



7-9 _ 1-75 ( 1 Tl — 1— ys ( I T'i — 1-75 , , 

^1 - — 1 I > -^^2 - — 1 5 ) , ^3 - — U I > (8) 



Tjt _ 1+75 „ 1+75 ,, r>l _ 1+75 _ 

r,q _ 1+75 „, n9 _ 1+75 „ pS _ 1+75 ^ 

-'^1,1 — 2 -"-1,2 — 2 -"-1,3 — 2 

-^2,1 ~ -^2,2 ~ -^2,3 ~ ^- 



The fermion part of the action is: 



Sf = J d^x {C\ + £f + C\ + £^'1 + + £Lss + -C^ssl- (9) 



Here, 



4 = iL\{d^ + iA^-iB^,)-f^Ll 

Cf^ = iRl^id^ + iZ^ + {4i/3)B^)j^Rl., 

= iRli{d^+iZ^-{2i/Z)B^)^^Rl., 

■^mass = lY.'^i{L'ir'^-R\. + \y.Mii{L'irn'^Ri.+h.c., 



V ■ 



mass 



i^mf(Zfr$«Ef + h.c.. (10) 

In these expressions = ia^^ (iui is the charge conjugation operator), ■i/i = V'Wo, and 

= me, m2 = m^, rrig = TOx, 

= m„, = mc, m| = mt- (11) 
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M is the mass matrix, which eigenvalues represent the masses of the d, s, and h quarks. The nondiagonahty of this 
matrix gives rise to the phenomenon of quark mixing. 

All necessary information about the Euclidian dynamics of the SM is contained in the gauge invariant correlators: 

(O(fields)) = j DZ DADBDeDeDveDve . . . L»$ exp(-S'(fields)) O(fields) (12) 

B. Representation of the Standard Model in loop space 

The hidden symmetry we are talking about may be seen after reformulation of the Standard Model through loop 
variables. (For the definition of the notations connected with loop space dynamics see Q.) The derivation is as 
follows. 

First, we note that in Eq. (|10|) L and R can be treated as independent two-component Weil spinors. In the Weil 
basis of the 7 - matrices the Euclidean fermion Lagrangian contains quadratic terms like L''^(Vo — i^i(Ti)L and 
-R^(Vo + iWiai)R (where V^^ is a covariant derivative, at {i = 1, 2, 3) are Pauli matrices), and an interaction term like 
(Lt$)i?. 

For an arbitrary gauge invariant correlator we have: 

(O(fields)) = J DZ DADB cxp{-Sg){0{^e\ds))f,s- (13) 

Here, (O(fields)) s is an integral over fermions and over the scalar field. It is calculated in the model with an external 
gauge field. First we perform an integration over the Grassmann variables. We can do it using simple Feynman rules. 
The diagrams contain propagators of Weil spinors, correlation functions of the scalar field, and interaction vertices 
that are coming from the term Lt"$"it;. The loops coming from the fermion determinant should also be taken into 
account. We use the path integral representation of the propagators and of the fermion determinant (see, for example, 
P, ) . In order to calculate the scalar field correlators we use the lattice regularization. It is well - known that the 
bosonic path integral for any field correlators in the external gauge field on the lattice has a representation as a 
sum over all possible closed loops. (For the details of the calculation see, for example, 0). After returning to the 
continuum representation we arrive at a path-integral representation of the scalar field correlators. 
Finally we represent any correlator Eq. H12I) in the form: 

(O(fields)) = J DCaO{C) J DZ DADB eM-Sg)W{C) =^ 

= f DCo.O{C){W{C)). (14) 



Here, Ca stands for the set of paths. Each path corresponds to one of the fermions (left - or right - handed) or to 
the scalar. It may either be closed or ending in a vertex. Each vertex corresponds to the transformation of left-handed 
fermion into right-handed ones and emission or absorption of a scalar. The definition of the measure DCa comes from 
the path integral representations of the bosonic correlator mentioned above, the fermion determinant, and the fermion 
propagator. It includes all possible paths described above. The index a enumerates all fermions and the scalar. The 
functionals O do not depend upon the gauge fields and are, hence, not of interest to us. The full dependence on 
the gauge fields is now concentrated in the loop variable W that is simply a product of parallel transporters W(Ca) 
corresponding to the fermions and to the scalar: 



W{C) = TyUWiCa). 



(15) 
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In Eq. (|15|l we encounter six different parallel transporters: 



/ 1/3 

Pexp(i / {Z^ + Af, + \B^,)dx^) = uj gu {3)^ su {2)^ u (^ly 
J c 

J c 

Pexp(i J^iZf, - ^Bf,)dXf,) = ujsui3)^^u^i) , 

Pexp(z y"(^^+ B^)dx^) = uJsu{2)i^u(i), 

where we introduced Wilson loops corresponding to SU{3), SU{2), and U{1) gauge fields respectively: 

^uii) = exp(i J Bf^dXfj,), 

'^SU{2) = Pexp{i J Af^dXf,), 

^su(3) = Pexp{i / Zf^dXf,). 

Jc 



WniiC) 
Wh{C) 



(16) 



(17) 



In Eq. H16|l each uj corresponds to a path connecting different points. However, in Eq. (|15() these parallel transporters 
are arranged in such a way that yV(C) depends only upon uj's corresponding to closed loops constructed of Cq. 

Thus any correlator can be represented through vacuum averages of products of those loop variables. The vacuum 
average is considered in the pure gauge theory with the action Sg. Using the loop calculus, we can express (yV(C)) as 
follows: 



(W(C)) 



Duju(i)Dujsu{2)Dujsu(3) x 



X exp - 



-Tr 



SU{3) 



5uj 



5(7(3) 



1 SuJsu{2) 



235(7(2) Sa^^ix) 

1 ^^U{1) 



Suj 



5(7(2) 



6a^^{x) 



<5crpy(x) 



5(Ti^,,{x) 



W(C), 



(18) 



where u^i, is the infinitesimal area, — is the area derivative, and . . . |o means that the area derivatives are calculated 
for infinitesimal contours. The measure over the gauge variables is denoted now as Duj. 



C. The symmetry 

Now we are in a position to point out the mentioned symmetry. It turns out that yV(C), being expressed through 
oj's corresponding to closed loops, is invariant under the following transformation: 

^u(i){C) exp(-i7rL(C,I]))tJ,7(i)(C), 
^5(7(2) (C) exp(j7rL(C,I]))cJs(7(2)(C), 

^5r/(3)(C) ^ exp(i|^L(C,E))w5t/(3)(C). (19) 



Here, S is an arbitrary closed surface and L(C, E) is the integer linking number of this surface and the closed contour 
C. From Eq. p9|) it is clear that this transformation belongs to Zq group. 
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This transformation corresponds to the centers of the SU{3) and SU{2) subgroups of the gauge group. It is finite, 
being apphed to the gauge invariant loop variables lo. However, it becomes singular in terms of gauge potentials: 



where Vfi{x) = Jytp_S{x — y{a,b,c))dadhdc is an integral over the three-dimensional hypersurface y(a,b,c), the 

boundary of which is S. The normal vector to V is denoted by tfj, = ^^fiupcr^g^^§^^§^- 

The invariance of Eq. p5|) under the transformations p9fl can be easily proven via direct substitution of Eq. 1)20(1 
into Eq. l(TB|l . 



D. Redefinition of tlie gauge action 



The whole SM can be represented (at least, formally) in such a way that it possesses the symmetry with respect to 
transformation Eq. ((19|l . This can be done by the following redefinition of the pure gauge part: 



(W(C)) 



J Duju{i)DuJsu(2)Di^su{3)ex.pi^Pk J Tr 

\ k 



5W^ 

{x) 



5Wk 




Here the sum is over the six parallel transporters mentioned above. For an appropriate choice of couplings Pk the 
action in Eq. (|21|l is equal to the action in Eq. (|18|l defined in terms of smooth gauge fields. However, in loop calculus 
we are not forced to consider smooth gauge fields. We are allowed to consider piecewise smooth loop variables u instead. 
The main difference is that the action in Eq. H18(l suppresses step-like w's, while Eq. (|21|l allows the appearance of 
loop variables with Eq. H19|) -like discontinuities. 

Instead, Eq. H21|) suppresses the discontinuities in Wk- Therefore we may apply an Eq. I|19|) -like transformation to 
all cj's in order to make them smooth. After that Eq. H21|) becomes identical to Eq. (|18|l . So, these two formulations 
would define the same theory. 

Here we implied that if the action suppresses some physical quantity, the latter indeed vanishes. However, there is 
another point of view. Namely, there are some indications, that the naively suppressed quantities may survive due 
to the entropy factor j^. We do not discuss here this possibility, but we must mention that if this picture emerges 
in the Standard Model, Eqs. (|18|l and (|21|) may define different models and correspond to different physics. We also 
notice here that in this case, say, the topological theta-term with 9 = 2tt being added to the action could, in principle, 
change the nonperturbative behavior of the theory while keeping the same perturbation expansion . 



E. The Standard Model as a finite cutoff theory 



It was mentioned in the introduction that the Standard Model should be regarded as a finite cutoff theory. So, the 
correct continuum model must contain a short distance part (related to the unification of the Electroweak and strong 
interactions), which makes the corresponding lattice model cutoff-independent. 

The Unified model could be the origin of our additional symmetry. If so, Eq. I|19|) emerges in it without any 
singular transformation of gauge potentials. Actually, in the corresponding examples considered in Q the realization 
of Eq. (|19|l being written in terms of the continuum fields is not singular. 

Strictly speaking, the only thing we are able to consider is the regularized model. The regularization can be 
constructed in such a way that it either admits or does not admit Eq. H19|l . Each choice of regularization is, in 
essence, the low energy limit of a regularized Unified model. After the discretization is removed, the full continuum 
theory appears. The finite cutoff Standard Model is an approximation to this hypothetical theory. Our assumption is 
that if we construct the finite cutoff Standard Model (FCSM) either keeping or not keeping the additional symmetry, 
the degree to which the resulting model approximates the correct unified model could be different. 

We expect that this difference might manifest itself at high enough energies. Probably, this could happen in 
the intermediate region between the usual SM scale and the GUT scale. Strictly speaking, in this region neither 
realizations of the FCSM can describe the physics properly. However, if the symmetry with respect to Eq. H19|l is 
indeed a fundamental symmetry, the corresponding model may give results that are closer to the experimental ones 



6 



(and vice versa). If so, we would catch the echo of the Unified model already at intermediate energies and draw 
certain conclusions about its structure. 

However, we expect that the most important role of the symmetry with respect to (|19|l is rather technical. The 
convergence of the lattice methods to physical results could become considerably faster for the models that respect 
invariance under Eq. H19() . This can be crucial for consideration of certain processes. Probably, the same situation 
takes place, say, for the SU{2) and 5*0(3) gauge models 0. They are generally beheved to belong to the same 
universality class. However, physical results are practically not achievable via 5*0(3) lattice theory. The reason is 
that the Z2 symmetry is lost. 

In general, it is thought that the convergence of a lattice model to the continuum results is faster if it keeps as much 
symmetries of the continuum model as possible. It even might occur that models that keep or do not keep a certain 
symmetry may lead to different continuum theories. Therefore, we also do not exclude that FCSM's that respect or 
do not respect Eq. H19|l would give essentially different results. In any case nothing definite could be said until the 
corresponding numerical research is performed. 



III. THE LATTICE MODEL 



A. Discretization of the continuum model 



In the remaining part of this article we shall not be interested in a discretization of the fcrmion sector. We would 
only notice that there are some difficulties concerning the problem of keeping the chiral symmetry while avoiding 
doubling. There were many different papers on this subject. For a review see |11.| and references therein. 

Now our aim is to remind the construction of 1]. We construct a lattice 5t7(3) x SU{2) x U{1) gauge model coupled 
to the scalar field in such a way that it reflects all the required properties of the Weinberg - Salam model and, in 
addition, preserves the symmetry considered above. 

The model contains the following variables: 

1. Lattice gauge fields (which live on the links of the lattice): 

r e 5C/(3), U e SU{2), e''^ e u{i). (21) 

2. A scalar doublet a = 1,2 (which lives on the lattice sites). The action of the model must have the form: 

S = Sg + Sh, (22) 

where we denote by Sg the pure gauge part and the scalar part of the action is denoted by Sh- 
A possible choice of Sh is 

Sh^YI lU.ye-'^^'^y - 3>.r + E (23) 

xy X 

where V{r) is the potential, which has a minimum at a nonzero value of r = ^/j. 

To construct the pure gauge part of the action we use the following correspondence between lattice and continuum 
notations: 

^S(7(3)(C) Hi (24) 

where I is a closed contour on the lattice corresponding to the continuum contour C. 
The analogue of the continuum transformation is the lattice transformation: 

U -> C/e"*'^^, 
e ^ 6l + 7riV, 

r -> re(2W3)w^ (25) 



where A'' is an arbitrary integer link variable. It represents a three-dimensional hypersurface on a dual lattice, the 
boundary of which corresponds to S in Eq. H19|l . This symmetry reveals the correspondence between the centers of 
the SU{2) and 5'?7(3) subgroups of the gauge group. 
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The choice /3 = = /Jr^ = Phg = /3ri = /3r2 and (3h = corresponds to a certain class of unified models Q. 
Hence, we choose: 

Sg = p J2 (2(l-^Tri7pCOS0p) + 

plaqucttcs 

+ (1 - cos26'p) + 

+6(1 - iReTrrpTrC/pexp(i6lp/3)) + 
+3(1 - iReTrrpexp(-2i6lp/3)) + 

+3(1 - iReTrrpexp(4i6'p/3))), (26) 

where the sum runs over the elementary plaquettes of the lattice. Each term of the action Eq. (|26|l corresponds to 
a parallel transporter along the boundary dp of plaquette p. The correspondent plaquette variables constructed of 
lattice gauge fields are Up = 0Jsu{2){dp) = i^su{3)idp) , and dp = ArgWf/(i)(9p). 

B. The simplified model 

In this paper we report our results on the numerical simulation of the model, in which we omit the dynamical 
fermions, as well as the color subgroup SU{3). It will be seen below that already on this level certain qualitative 
differences between this model and the conventional one exist. 

The potential for the scalar field is considered in the London limit, i.e., in the limit of infinite bare Higgs mass. The 
action of the model reduces to 

iS* = + + Sh ~ 

= P ((l-iTrf/pCos6»p) + i(l-cos26'p)) + 

plaquettes 

+ Y,\U.ye-'''''%-<fx\^ + Vm). (27) 

xy 

(Here f3 is rescaled as /? ^ f3/2 for the convenience of comparing the results with those of the SU{2) fundamental 
Higgs model.) corresponds to the doublet of left-handed fermions and corresponds to the right-handed singlet. 
$ is the Higgs doublet and V is an infinitely deep potential, giving rise to the vacuum average (|$|) = y^. It is 
worth mentioning that the naive continuum limit of Eq. I|27() gives the value of Weinberg angle 0w = 7r/6, which is 
surprisingly close to the experimental value. 
After fixing the unitary gauge we obtain: 

S = 13 ((1 - i Tr [/p cos 6'p) + 1(1- cos 26'p) + 

plaquettes 

+7E(l-^e(C^iie-'^'')))- (28) 

xy 

Of course, we keep in mind that this simplification of the model may lead to some qualitative changes in the description 
of the dynamics. Thus the conclusions, which we draw after performing the numerical investigation of the simplified 
model, must be justified by the study of the full model, including the color subgroup, dynamical fermions, and a finite 
Higgs mass. 

Below we briefly describe some of the quantities which we investigate in this work. 

The following variables are considered as creating a photon, Z boson, and W boson respectively: 

A^y = A>i = [ArgC/ii+0,,]mod27r, 
Z^y = Z^ = [Argf/ii-0,j,]mod2^, 

W^y - Wi^ = Ulle^'^y. (29) 
Here, p. represents the direction {xy). After fixing the unitary gauge the electromagnetic U{1) symmetry remains: 

Uxy dxUxydyi 

Oxy e^y + ay/2-a^/2, (30) 
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where gx = diag(e*"'"/^, e '"^^Z^). The fields A, Z, and W transform as follows: 



A 


^ -^xy ^y ^xt 


y 

^xy 


^xyi 


Wxy - 





(31) 

As any other compact gauge theory, our model contains monopoles. As in other compact gauge models, their 
behavior is connected with the possible confinement of charges. On the other hand, the continuum Weinberg Salam 
model is believed not to confine any charges and not to be affected by monopoles. 

We investigated two types of monopoles. J7(l) monopoles extracted from 20 are defined as 

ne = 7^*d([d26']mod27r). (32) 
27r 



The electromagnetic monopoles are: 



jA = -^*d([dA]mod27r). (33) 

ZTT 



(He re we used the notations of differential forms on the lattice. For the definition of those notations see, for example. 
The density of the monopoles is defined as follows: 

/ X] links blink I \ 

P \ 4X4 /' ^-^^^ 

where L is the lattice size. To understand the dynamics of external charged particles, we consider the Wilson loops 
defined in the representations of left-handed and right-handed leptons: 

>V^(0 = (ReTrn(,j,)e,[/,j,e-'^-), 

WJ^(0 = (Ren(,^)e,e-2^«-). (35) 

Here I denotes a closed contour on the lattice. We consider the following quantity constructed from the rectangular 
Wilson loop of size a x a: 

Vi?,L (a) = - log W^'^ (a xa)/a. (36) 
A linear behavior of Via) would indicate the existence of a charge - anti charge string with nonzero tension. 



C. Numerical results 



In our calculations we investigated lattices for L — 6, L — 12, and L — 16 with symmetric boundary conditions. 

We summarize our qualitative results in the phase diagram represented in Fig.^ The model contains three phases. 
The first one (I) is a confinement-like phase, in which the dynamics of external charged particles is similar to that of 
QCD with dynamical fermions. In the second phase (II) only the behavior of left-handed particles is confinement-like, 
while for right-handed ones it is not. The last one (III) is the Higgs phase, in which no confining forces are observed 
at all. This is illustrated by Figs. [3 and El in which we represent VL{a) and VR{a) at three typical points that belong 
to different phases of the model. One can see that in the Higgs phase the shape of V(a) excludes the possibility of 
a linear potential to exist. The same behavior is found in phase II for V_R(a). On the other hand, in phase II the 
shape of VL(a) signals the appearance of a linear potential at sufficiently small distances (up to five lattice units). 
However, as for QCD with dynamical fermions or the SU{2) fundamental Higgs model 13, ^lij, these results do not 
mean that confinement occurs. The charge - anti charge string must be torn by virtual charged scalar particles, which 
are present in the vacuum due to the Higgs field. Thus V(a) may be linear only at sufficiently small distances, while 
starting from some distance it must not increase, indicating the breaking of the string. Unfortunately the accuracy of 
our measurements does not allow us to observe this phenomenon in detail. However, it may be partially illustrated 
by the shapes of Vi(a) and Vi{(a) in phase I shown in Fig. Inland Fig. 13 

The phase structure of the model may also be seen through the data for the mean action over the whole lattice 
S — {S) /{6PL^), Fig. ^ It appears to be inhomogeneous in a small vicinity of the phase transition line. 
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The connection between the properties of monopoles and the phase structure of the model is illustrated by Figs. [S] 
and El which show the monopole density versus the coupling constants. The electromagnetic monopole density drops 
in the Higgs phase, while the t/(l) monopole density falls sharply both in phase II and phase III. We can see that the 
behavior of the U{1) monopoles is connected with the dynamics of the right-handed particles, while the behavior of 
electromagnetic monopoles reflects the dynamics of the left-handed particles. 

It is worth mentioning that the cousin of our model, the SU (2) fundamental Higgs model, has a similar phase 
structure as our model, except for the absence of the phase transition line between phases I and II. In the latter 
model it was shown that different phases are actually not different. This means that the phase transition line ends 
at some point and the transition between two states of the model becomes continuous. Thus one may expect that in 
our model the phase transition line between phases I and III ends at some point. However, we do not observe this for 
the considered values of couplings. 

In our model both phase transition lines join in a triple point, forming the common line. This is, evidently, the 
consequence of the mentioned additional symmetry that relates SU{2) and U(l) excitations. The same picture, of 
course, does not emerge in the conventional SU{2) (g) [/(I) gauge - Higgs model 

We must also notice here that the phase diagram may also contain an unphysical region, corresponding to the 
unphysical region of the pure SU(2) model (which is observed at /? < /3c, where Pc is the crossover point). Our 
investigation shows that if this region of couplings exists in our model, it must be far from the Higgs phase, which is 
of main interest to us. Indeed, this unphysical region might appear for /3 < 2.25 and 7 < 0.5. 



IV. CONCLUSIONS 



We summarize our results as follows: 

1. We illustrated an additional symmetry found in the fermion and the Higgs sectors of the Standard Model by the 
consideration of the SM in loop space. 

2. We performed a numerical investigation of the quenched Electroweak sector of the lattice model, that respects 
the additional symmetry. 

3. The lattice model contains three phases. The first one is a confinement-like phase. In the second phase the 
confining forces are observed, at sufficiently small distances, only between the left-handed particles. The last one is 
the Higgs phase. 

4. The main consequence of the emergence of the additional symmetry is that the phase transition lines corre- 
sponding to the SU{2) and U{1) degrees of freedom join in a triple point forming the common line. This reflects the 
fact that the SU(2) and U{1) excitations are related due to the mentioned symmetry. The same situation does not 
take place in the conventional SU{2) (g) U{1) gauge - Higgs model [l5l| . 

So, already on this simplified level, we found a qualitative difference between the conventional discretization and 
the discretization that respects the invariance under Eq. (fl^ . 
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FIG. 2: Vi(a) calculated at three points that belong to different phases of the model. 
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FIG. 3: Vh(o) calculated at three points that belong to different phases of the model. 




FIG. 5: The density of electromagnetic monopoles. It decreases in the Higgs phase. 
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FIG. 6: The density of (7(1) - monopoles. It decreases when the behavior of right - handed external particles is not confinement 
- like. 



